Theory of Umklapp-assisted recombination of bound excitons in Si:P 



in 
o 
o 

(N 



in 



Michael N. Leuenberger and L. J. Sham 
Department of Physics, University of California San Diego, La Jolla, CA 92093 

(Dated: February 2, 2008) 

We present the calculations for the oscillator strength of the recombination of excitons bound to 
phosphorous donors in silicon. We show that the direct recombination of the bound exciton cannot 
account for the experimentally measured oscillator strength of the no-phonon line. Instead, the 
recombination process is assisted by an umklapp process of the donor electron state. We make use 
of the empirical pseudopotential method to evaluate the Umklapp-assisted recombination matrix 
element in second-order perturbation theory. Our result is in excellent agreement with the experi- 
ment. We also present two methods to improve the optical resolution of the optical detection of the 
spin state of a single nucleus in silicon. 

PACS numbers: 71.35.-y, 78.67.-n, 03.67.Lx, 76.70.Hb 
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I. INTRODUCTION 

Proposals on quantum computing in semiconductors 
have recently attracted a great deal of attention. ^'^ The 
main idea is to use the electron spin of quantum dots in 
semiconductors^ or the nuclear spin of shallow donors in 
silicon^i as a qubit for quantum information processing. A 
complete quantum computation consists of, besides the 
single- and two-qubit operations, the initialization and 
the readout of the qubits. While the initialization and the 
readout of the electron spin in quantum dots have been 
successfully demonstrated experimentally^£ it remains 
an experimental challenge to read out the nuclear spin of 
a donor electron in silicon. 

It has recently been proposed that the photolumines- 
cence of excitons bound to phosphorous donors can be 
used to detect the spin state of a single donor nucleus in 
silicon.^ This scheme for optical readout could render the 
recent quantum computing proposal using conditional 
NMR and ESR pulses feasible .-SJ^ Although experiments 
have shown that the recombination of the bound exci- 
ton follows strong optical selection rulesjii it has been 
unclear what physical process is responsible for the no- 
phonon line of the optical recombination of the bound 
exciton. In contrast to the phonon-assisted recombina- 
tion, the no-phonon line represents the recombination 
process without phonon assistance. A shell model that 
accounts for the selection rules has been proposed in 
Ref. fl^ T he shell model has been later improved in 
Ref. Il3ll4l by a Hartree-Fock calculation that takes the 
multivalley character of indirect bandgap semiconductors 
into account and therefore is in good agreement with the 
measured fine-structure excitation spectrum of the bound 
exciton complexJ^. The Hartree-Fock calculation fails to 
predict binding energies. This problem was overcome 
by a density-functional calculation that takes the cor- 
relation energy into account.^^ However, neither model 
gives a satisfactory physical description of the recombi- 
nation process for the no-phonon line and thus is un- 
able to quantitatively reproduce the measured oscillator 
strength /cxp in Ref. |0 Our calculations show that the 
probability for direct recombination of the bound exci- 



ton is negligibly small. Here we present a physical model 
of the recombination process that accurately reproduces 
the oscillator strength /oxp of the no-phonon line. In our 
model the recombination of the exciton is assisted by the 
Umklapp-process of the donor electron. In Sec. |n] we 
give a detailed description of our model, where we intro- 
duce the Coulomb scattering of the bound and the donor 
electron. For the evaluation of the recombination matrix 
element, we make use of the empirical pseudopotential 
methodi24i2i2£ to calculate the bandstructure of silicon 
with 137 reciprocal lattice vectors at 100 points inside 
the first Brillouin zone in X direction. This technique 
is applied to our model in Sec. Hn Then we use the re- 
sulting Bloch states to calculate the oscillator strength of 
the Umklapp-assisted recombination in second-order per- 
turbation theory, which is shown in Sec. IIVI In Sec. El 
we present two methods based on optically detected mag- 
netic resonance (ODMR) to improve the resolution of the 
readout of a nuclear spin of a donor electron in silicon. 



II. PHYSICAL MODEL 

Our physical model of the Umklapp-assisted recombi- 
nation is shown in Fig. ^ First the bound electron with 
wavevector ke is scattered via the Coulomb potential 
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off the donor electron with wavevector Icd, thereby 
conserving the momenta, i.e. ke + kn = kj 
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eo = 8.854 x 10~^^ F/m is the dielectric constant of 
the vacuum, and = 11.56 is the relative dielectric 
constant of silicon, which are combined to e = eoCr^ 
So the refractive index of silicon is about i> = y/e^ = 
3.4. £,1 = Amikpq^/'nh^ = 1.9 x 10^ m^^ and = 
imtkFq^ /nfP = 4.0 x 10*^ m^^ arc the Thomas-Fermi 
screening lengths in longitudinal and transverse direc- 
tion, where kp = (STr^rto)^^'^ is the Fermi wavevector 
and no = is the density of the donor electronsi^ 

The distance between the donor electrons is about L = A 



2 



nmm^ The longitudinal mass in silicon is mi = 0.9163me, 
where rrie = 9.1095 x 10"^^ kg is the bare electron mass. 
The transverse mass in silicon is mt = 0.1905me. After 
the Umklapp process the bound electron with wavevec- 
tor kg rccombines with the bound hole with wavevector 
kj^ = kg via the electron-photon interaction 



mt 



(2) 



which results in the emission of a photon that can be 
seen in a photoluminescence experiment. So the total 
Hamiltonian describing our system is given by 



n 



^ 2m* 



+ Vc + V„ 



(3) 



in a spin singlet state (ITbId) - |iBTD))/V2, the or- 
bital wavefunction is symmetric in the valley combina- 
tions j = ±ka;o,±kyo,±k^o, i.e., 



BD 
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(5) 



(0,±A:o,0), ±k,o 



1 _ H + h 



where ±'kx:o = (±/co,0, 0), ik^o 
(O,0,±fco), and e.g., 

Fek.o(r)-J^Bk.o(r) = i^Dk.o(r) 

(6) 

is the envelope function of the bound and the donor elec- 
tron with a = 25.1 A and b = 14.4 A (see Refs. ,24^2h) . 
and 



where the first term is the kinetic energy of the bound 
electron, the donor electron, and the bound hole with the 
effective masses m* , which takes the lattice potential via 
the empirical pseudopotential method into account (see 
App.m). 



0ok,o(i") = 0Bk,o(r) = '/'Dk,o(r) = u^z{Tc)e'''°' 



(7) 
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FIG. 1: Although the Coulomb scattering amplitude Aic has 
a maximum at ku = 0.55fcsi (see Fig. I^J, which arises from 
the scattering process (b), the largest contribution to the 
Umklapp-assisted recombination comes from the scattering 
process (a), because the overlap of the unbound scattered 
electron and the hole is largest around ku ~ (see Fig. 0)1. 

The transition amplitude of the Umklapp-assisted re- 
combination is given by 



1 1 
Aftot = T / dkuM{k^) = — 

ko Jo Ko 



dk„ 



McMr 



inor -^c(^u) 

(4) 

in second-order perturbation theory. -Edonor is the en- 
ergy level of the donor electron and i?c(fcgo) is the en- 
ergy of the conduction band. In order to calculate the 
Umklapp-assisted scattering matrix element Mq and the 
recombination matrix element Mr, we need to identify 
the initial, intermediate, and final states (see Fig. ^ and 
App. 0. Since the donor and the bound electron are 



is the Bloch wavefunction. ko — 0.85fcsi is the distance 
from the F point to the minimum of the conduction band 
in X direction (see next section), kgi — 2iT/d is the re- 
ciprocal lattice vector of silicon with a lattice constant of 
d = 5.43 A. Similarly the envelope function of the hole 
state is given by 



J^h(r) 



1 



(8) 



with c — Aireh? /miiY^q^ , where mhh = 0.523me. The 
Bloch wavefunction of the hole is 



<^h(r) = Mv(r)- 



(9) 




FIG. 2: The Coulomb scattering amplitude between the 
bound and the donor electron has a maximum at ku — 0.55/Csi. 

For the Coulomb scattering matrix element we obtain 
Mc = '^iC-^ / dkBz / dk^z^'^oif^Bz - K + kuz) 



1 



{kBz ~ ku) + 



-i'Bo{kBz)'4'Uoik-Dz), (10) 



which is shown in Fig. |21 The maximum of Mq is located 
at ku = 0.55A:si, which is due to the scattering process (b) 
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shown in Fig. ^ For the recombination matrix element 
we obtain 



hqAi 



<fKi:^\K)KMK). (11) 



which is shown in Fig.^ M-^ has a maximum at fc„ = 0. 
The analytical derivations of Mq and Mr can be found 
m App.EI Next we need to calculate the bandstructure 
of silicon, which is presented in the next section. 



III. BANDSTRUCTURE OF SILICON 
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FIG. 3: Bandstructure of silicon in the direction of the X 
point. The minimum of the conduction band with Ai sym- 
metry is located at fez = 0.85A;si. 

We choose the empirical pseudopotential method to 
calculate the bandstructure of siliconii2ii2iSi A brief re- 
view of this method is given in App. El We need to 
solve the Schrodinger equation for the pseudowavefunc- 
tion, i.e., 



n^(k + G)^ 
2me 



G' 



)UG'=EUG. (12) 



We diagonalize it using the 137 G-vectors shown in Tab.^ 
and the form factors Vs{G = V3) = -3.049 eV, Vs{G = 
V8) = 0.750 eV, and Vs{G = 711) = 0.985 eV. 

We calculate the bandstructure of silicon for 100 points 
in k-space in the direction of the X point. The band- 
structure is shown in Fig. |31 For our Umklapp-assisted 



TABLE I: G-vectors used in the empirical pseudopotential 
calculation. 



recombination the following symmetries are important: 
the conduction (valence) band at the F point has the 
symmetry r25' (Fis), which transforms as {yz,xz,xy} 
{{x,y, z\). The conduction band at the X point has the 
symmetry Ai, which transforms as {z}. So the initial 
bound electron state before the recombination has the 
symmetry Ai and the hole has the symmetry Fis. Since 
the spin-orbit splitting at the F point is Ago = 0.044 eV 
(see Ref. |2^ . the heavy-hole states including spin are 
represented by 



J = -:M : 

2' 



■^{\yz)+i\xz))\]), 

2 1 

3 '^^ ~ 71 \{y + ix)z i) , 

2 1 

3 ^'^^vE 1(2^" ' 

i=(|yz)-zM)|i). (13) 



That is the origin of the selection rules found experimen- 
tally in Ref. 11. 
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FIG. 4: The recombination amplitude for the unbound scat- 
tered electron is largest around ku = 0. The oscillator 
strength for the recombination of the scattered electron is 
/rcc(feu = 0) = 0.023. 
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IV. OSCILLATOR STRENGTH 

In textbooks such as Ref. QJi the oscillator strength is 
defined as 



mEn 



\{n\p\i)\\ 



(14) 



where Eni is the energy difference between the initial 
\i) and the final state |n). For our Umklapp-assisted 
recombination the oscillator strength is thus given by 



^opt -^gap 



ghhQsMtot 



rrit 



(15) 



where -Egap ~ 1-1 cV is the indirect bandgap energy, 
iTiopt — 3 / {1 / mi + 2 / nit) is the effective mass of the elec- 
tron that ensures the sum rule J2n fm = Ir^ ghh — 2 
is the degeneracy of the heavy-hole band and (?s = 2 is 
the spin degeneracy. We calculate first numerically the 
Coulomb scattering amplitude Mc{ku), which is shown 
in Fig. 121 Then we calculate numerically the recom- 
bination amplitude Mr{ku), which is shown in Fig. 0] 
In Fig. [3 we plot the Umklapp-assisted recombination 
amplitude M{ku)- Inserting Mc(fc„) and Mriku) into 
Eq. (0} and integrating over ku yields 



/tot = 10.1 X 10- 



(16) 



which is in excellent agreement with the oscillator 
strength /oxp = 7.1 x 10^^ of the recombination of an 
exciton bound to a phosphorous donor in silicon reported 
in Ref. llTl For comparison, the direct recombination of 



the bound electron with the bound hole over the indirect 
bandgap has an oscillator strength of /direct = 4 x lO^'^'^. 



V. OPTICAL READOUT 

In this section we show two methods that can be used 
to improve the optical resolution of the optical detection 
of the spin state of a single nucleus in silicon (see Ref.0): 

• The first method is inspired by the optically de- 
tected magnetic resonance technique (ODMR)«, 
where the induced ESR Rabi oscillation alters the 
lifetime of the bound exciton and is thus optically 
detectable. Since the transition from the bound 
hole state Mj = -3/2 to Ms = +1/2 is forbidden, 
mixing the spin states in an equal superposition of 
Ms = -f 1/2 and Ms = -1/2 leads to a doubling of 
the lifetime of the bound exciton. 

• The second method makes use of a strong ESR field 
that renormalizes the spin levels in such a way that 
the hyperfine and Zeeman splitting of the nuclear 
spin is increased. As in the first method, the pho- 
tons of the ESR field dress the spin states of the 
donor electron. 

The Hamiltonian for the electron-nucleus system is 



Hs = geMsH • S - ,g„^„H ■ I + AS • I 



(17) 



in the generalized rotating frame. If we choose the basis 
{iTetn) , |Tei«> , UeT«) , |iei«)}, the Hamiltouian reads 
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The energy levels of the electron and nuclear spin are We follow here the derivation of dressed states given in 
shown in Fig. El Ref. The Hamiltonian of the quantized photon field 

is 



Typical external magnetic fields used in experiments 
are i?^ = 3.0, 5.0, and 7.0 T. In a field of H, = 3.0 T 
the Zeeman splittings of the electron and nuclear spin 
are ge^J^BHz = 0.35 meV and gnl^nHz — 0-21 /ieV, re- 
spectively. The hyperfine splitting is A — 0.50 /leV. 
Thus the splitting between |teTn) and Itein) is Af^ = 
-^TeTn ^ ^Tei,i = ^/2 — gnf^nHz = 0.03 /^eV, whereas 
the splitting between lietn) and |iein) is = ^^iei„ - 
-^ieT" — + .9riMn-ffz — 0-46 fJ.eV. The hyperfine and 
Zeeman splittings for — 3.0, 5.0, and 7.0 T are shown 
in Tab. ini 



Hp = hua^a, (19) 

where d (a^) is the annihilation (creation) operator of 
a photon. The interaction between the microwave pho- 
tons and spin of the electron of the Phosphorus donor is 
described by 

V = -n B, (20) 
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FIG. 5: Due to the large overlap of the unbound scat- 
tered electron and the hole around fc„ = (see Fig. 21 1 
the Umklapp-assisted recombination amplitude Mtot is dom- 
inated by the terms M{ku) around ku ~ 0. For example, 
/(fc„ = 0) = 1.1 X 10"^ 



where the quantized magnetic field is given by 



B 



i(kx e) ^. (-i) (k x e*) 
k 



, (21) 



where Cp is the hght velocity. The spin-photon interaction 
can be simplified to 



V = gsp [(e • S) a + (e* • S) at] 
We use the circular polarization vectors 

e± (e^ ± iey) 



(22) 



(23) 



and the spin ladder operators S± = iSy. Then we 
obtain 



v2 



There are two coupled states 



= \ie,N+l), 
l0Te> = ITe,A^). 



(24) 



(25) 



If the interaction vanishes, the energies are E^^ = {N + 
l)hLLj — ^hujQ and = Nhui + ^hajQ. The energy sep- 
aration is huji^i^ — h{uj ~ luq). The matrix elements of 
V„+ read 

{4>u\v,^ = (0TelK.+ I0tJ = 0, 

The eigenstates are 



\Xi{N)) = sin0|<^^J+cos0|0tJ, 

|X2(^)) = cos0|0xJ-sin0|0tJ, (27) 



where tan 261 
ergies are 



El 



/2 
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Q./ (uj - ujq), <20 < it. The eigenen- 

(28) 



hid ± h\ 



LUq 



Let us tune the oscillating transverse microwave field to 
the transition between |teTn) and lietn)- Then the eigen- 
states are 



\XU2) 

where tan26' 
electron-nucleus system are 

9etJ-B 



(l0i.) + l0T.))ITn)/%/2, 

ism9\<PiJ+ COS e\cj)^ J) \U)/V2, 
{cos0\(t>iJ-sm0\^^^))\U)/V2, 
i\<t>u)-\<f>u))nn)/V2, (29) 

— ^gef^sHx/A. The eigenenergies of the 
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2 4 2 



E\,a 



The oscillator strength for the exciton recombination in 
Si:P is / = 7.1 x 10"^ (see Ref.[l3). The binding energy 
of the exciton to the phosophorus donor is i?binding = 4.7 
meV. The recombination rate is about w — 400 s^^. This 
leads to an interaction energy of 




E;nt = Vfiwf = 12 neV, (31) 
where V = 100 /xeV is the linewidth of the recombination. 







Qn f-^n Hz 






3.0 T 
5.0 T 
7.0 T 


0.35 meV 
0.58 meV 
0.81 meV 


0.21 ^ieV 
0.36 ^leV 
0.50 ^leV 


0.03 

-0.11 /xeV 
-0.25 /ieV 


0.46 ^leV 
0.61 ^eV 
0.75 



TABLE 11: The hyperfine and Zeeman splitting of the spin 
states of the donor electron. 

If the electron-photon interaction energy i?int is much 
weaker than the microwave coupling energy g^figH^, the 
states |xT„i), \xi„i), IXi„2), and |xt„2) are good eigen- 
states. Then the ODMR method works, because the final 
state of the radiative recombination is given by |xTni)j 
which doubles the lifetime of the bound exciton. The 
transverse microwave field should have a strength of at 
least Hx = 1.0 G, leading to a microwave coupling energy 
of geiisHx = 10.0 neV. 

In addition, it can be seen from Eq. H30|) that the 
hyperfine splittings are increased, which effectively in- 
creases the optical resolution. In order to obtain at 
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FIG. 6: The energy levels of the electron and nuclear spin of 
the donor electron are shown at the botton. The energy levels 
of the bound hole are shown at the top. 



least a 10% increase of the hyperfine+Zeeinan splitting 
of the two lowest energy levels, the oscillating trans- 
verse magnetic field must have a strength of 10.0 G, 
which yields the transverse Zeeman splittings of the 
electron and nuclear spin of ge^Jt-sHx — 0.12 /leV and 
gnUriHx = 0.07 neV, respectively. For a longitudinal 
magnetic field of = 3.0 T, the shifted hyperfine 



E 



Tn2 



E 



U2 



0.10 /ieV, whereas 



splitting is A2 
Ai = Ei^^i - Ef^^i = 0.53 ^leY. The shifted hyperfine 
and Zeeman splittings for Hz — 3.0, 5.0, and 7.0 T are 
shown in Tab. IIIII 

The frequency-resolved photon detection has the ad- 
vantage that a single photon detection is sufficient for 
determining the spin of the nucleus, whereas the ODMR 
measurement needs to be done in an ensemble in order to 
determine the radiative recombination time. Maybe the 
hyperfine-|-Zeeman splitting of A|_. = A/2 + gnfinHz = 
0.75 ^eY in a Hz = 7.0 T field is already sufficient for 
frequency-resolved photon detection. 



VI. CONCLUSION 

We calculated the oscillator strength of recombination 
of an exciton bound to a donor electron in silicon. We 









A2 


Ai 


3.0 T 
5.0 T 
7.0 T 


0.03 ^eV 
-0.11 fieY 
-0.25 fieV 


0.46 AieV 
0.61 fieY 
0.75 fieV 


0.10 AteV 
-0.04 AteV 
-0.18 AteV 


0.53 A^eV 
0.68 A^eV 
0.82 AteV 



TABLE III: Change of the hyperfine and Zeeman splitting of 
the spin states of the donor electron. 



showed that the Umklapp-assisted recombination, con- 
sisting of a scattering between the bound electron and 
the donor electron and the recombination of the un- 
bound scattered electron, gives the main contribution. 
The calculation of the Umklapp-assisted recombination 
was done in second-order perturbation theory. We made 
use of the empirical pseudopotcntial method to find the 
Bloch wavefunctions of silicon. The calculated oscillator 
strength is in excellent agreement with the experiment. 
We gave also two methods to improve the resolution of 
the detection of a nuclear spin of a donor electron in sil- 
icon. 
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APPENDIX A: MATRIX ELEMENT OF 
UMKLAPP-ASSISTED RECOMBINATION 

We are going to calculate the matrix elements Mc 
M-c in the reciprocal lattice space by Fouriertransform- 
ing the wavefunctions. The Bloch wavefunctions can be 
expanded in reciprocal lattice vectors as 



G 

for the electron wavefunctions and 



iGr 



(AI) 



(A2) 



for the hole wavefunction. Thus the electron 
wavefunctions in k-space are given by '0e(k, fcc) = 
'^cGV'eG(k, kc) with 



1 



1 



1 

+ {kz + Gz~ kef) 



(A3) 
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centered at fee, and the hole wavefunction by iphO^) 
I]G"^GV'hG(k) with 



1 



(^ + (fc, + G,)2) 

1 
1 



(A4) 



So the initial bound and donor electron wavefunctions are 
V'B(k) = V'D(k) = '!/'c(k, fco). The intermediate scattered 
electron wavefunction in k-space is given by ■!/'B(k') — 

I]G'"cGV'BG(k') with 



V^BG(k') = 



1 



1 



Gy 



(A5) 



where = lOOOd, the intermediate donor electron wave- 
function by V'i3(k') = Vc(k, /coo), where /cqo = 1.15/csi, 
and the intermediate hole wavefunction by i/j^k') — 
■0h(k). Note that only the bound and the donor elec- 
tron scatter off each other. 

Since the Coulomb interaction is local within each Bril- 
louin zone, interference effects can be neglected and thus 
it is sufficient to calculate the Coulomb scattering matrix 
element within the first Brillouin zone, i.e. 

Mc = {iP'^^'j.IVcIMb) J d^'k'^ Jd^kB j d^k^ 
xV''Bo(kB)V''Do(kB - k^ + ki 



f 



(kB-kJ,r+C2 
xViBo(kB)V'Do(kD), (A6) 

where we used the relation 



47rerBD 



(2^)3e 

„2 



3 (5(kB - k^ -I- kBD)'5(kD - k^ - kBD) 

^ ''BD p I ^2 

-<5(kB-kj3 + kD-ki3), (A7) 



(2^)3e(kB-kB)'+e- 



where tbd = tb ^ td and kBo = 27r/rBD- Numerical 
calculations show that the Coulomb potential is well ap- 
proximated by 



Vc(kB-ki3) 



f 



e {kB.-k'^,Y+e 



(A8) 



i.e. only the longitudinal dependence of the Coulomb 
potential in z direction is taken into account. Then the 
integrations in x and y directions can be solved analyti- 
cally, which yields 



Mc 



r'^ I dk' 



dk 



Bz 



X^'mif^Bz - fcL + ^Bz 
X^BoikBz)'ipI)o{kBz), 



dkBz^'*Bo{k'Bz) 



[kBz-kU +^ 



where 



1 



Tj_ = 



47r2 

X 



- U [log(- 
a I 



■la) 



log(ia) 



(A9) 



log(-j^) - log(i^) 



log(-i|) - log(*|) + log(-m) - log(ia) 
+i [log(-m) - log(m)]^ log(-j|) - log(i^) 



-71og(-za) + 71og(za)]} 



(AlO) 



Since the scattered electron state ■0'bo(^Bz) ™" 
bound, we can treat it as the square root of a 
delta-function with width i.e. ■'/''bo(^Bz) 



(5(<^f ) (fc^^ — ky). This means that we can approximate 

the integral / dk'Q^'4>'*Ba{k'^^)h{k'^^) by ft.(fc„)C, where 

C = J dk'^^tp'mi^Bz) ^ function h^k'^^). Thus we 

obtain 



Mc = tIC^ J dkBz 

1 



[kBz-ku) +C 



dkBz^'*Bo{kBz -ku + koz) 
i'Bo{kBz)'ipuoikBz), (Afl) 



which is solved numerically (see Fig. EJ. 

In the case of the recombination we cannot neglect 
interference effects, because the electron-photon interac- 
tion is nonlocal in k-space. We obtain 



(A12) 



which is solved numerically (see Fig. 



APPENDIX B: EMPIRICAL 
PSEUDOPOTENTIAL METHOD 

We give here a brief review of the empirical pseudopo- 
tential method. At each lattice site there is an atom with 
a nucleus, core electrons, and valence electrons. The at- 
tractive nuclear potential Vn is large and varies strongly 
throughout the lattice. The main observation is that Vn 
is almost canceled by the repulsive potential Vrcp pro- 
duced by the core electrons. So we need to consider only 
the valence electrons moving in a net weak one-electron 
potential VJ,, which is called the pseudopotential. 
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The full Bloch wavefunctions can be expressed as 

t 

$ must be orthogonal to the core wave functions (pt, i.e. 
((pt|$) = 0, which yields 



(B2) 



Then applying the Hamiltonian Hn = P^/'^rric + to <i> 
leads to the Schrodinger equation 



P 



2m, 



+ K + Vr 



rep 



E<t>, 



(B3) 



where the short-range non-Hermitian repulsive potential 
is given by 



(B4) 



with E being the full energy eigenvalue of $. ip is called 
the pseudowavefunction. Since the crystal potential is 
periodic, the pseudopotential Vp = Vn + Kcp and can be 
expanded in a Fourier series over the reciprocal lattice 
vectors G, i.e. 



G 



G 



where 



(B5) 



(B6) 



For zincblende and diamond lattices usually a two-atom 
basis is chosen, such that 

^0 (r ) = Vcation (r - r) + 14„ion (r + r) , (B7) 

where t = (1, 1, l)d/8. Then the Fourier potential reads 

Vb(G) = e''°-^Kation(G) + e-'°-^Knio„(G). (B8) 

The Fourier coefficients can be rewritten in terms of the 
symmetric and antisymmetric form factors Vs = Vacation + 

Kinion and Va = Vacation - Kinion. ThuS 

VoiG) = cos(G • T)Vs{G)+ism{G ■ t)Va{G), (B9) 

where the prefactors of the form factors are the structure 
factors. Since silicon has a diamond lattice, Va{G) = 0. 

The pseudowavefunction can also be expanded in a 
Fourier series, i.e. 



(r) = e'^'^u 



(BIO) 



G 



Inserting the pseudowavefunction (p and the pseudopo- 
tential Vp into Eq. IjBSp yields 



^-i^±^UG + y2Vo{\G-G'\)uG^=EuG. (Bll) 



G' 



Diagonalization of this Hamiltonian yields an effective 
mass Hamiltonian, which is the first term on the right 
side of Eq. ©. 
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